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the liner's ability to compress arising from the energy sink introduced by ionization 
and electron excitation, which are not accounted for in a polytropic EOS model. 
Variation of the plasma liner species for the same initial liner geometry, mass density, 
and velocity was also explored using the LTE tabular EOS model, showing that the 
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I. INTRODUCTION 



The Plasma Liner Experiment (PLX) at Los Alamos National Laboratory was designed 
to form and study spherically imploding plasma liners via the merging of thirty supersonic 
plasma jets.- Such imploding plasma liners could potentially allow repetitive generation of 
cm-, fis-, and Mbar-scale plasmas for fundamental high energy density plasma physics stud- 
ies, or (if scaled up in energy) be used as a standoff compression driver-"- for magneto-inertial 
fusion (MIF).-"- Recent one-dimensional (ID) radiation-hydrodynamic simulations of tar- 
getless imploding spherical argon plasma liners by Awe et ali^ used a polytropic equation-of- 
state (EOS) model, providing insight into the ID physical evolution of imploding spherical 
plasma liners, and showing that PLX-relevant plasma liner kinetic energies of hundreds of 
kJ could result in liner stagnation pressures on the order of 1 Mbar sustained for order 
1 /is. The results of Awe et ali^ indicated high enough plasma temperatures (lO's of eV) for 
ionization effects, which are not included in a polytropic EOS model, to have an appreciable 
impact on both the predicted peak pressures and temperatures. As some of the initial liner 
kinetic energy would now be used for ionization, it is expected that the liner may compress 
less, although the expected drop in thermal energy in the imploding plasma liner due to 
ionization and radiative loss channels are expected to help compensate. Note that this work 
and that of Awe et ali^ focus on targetless, self-collapse of the liner, motivated by the goals 
of the PLX project. Thus, care must be taken in using these results to provide insight 
into plasma liner implosions onto a magnetized target for MIF, which has been explored in 
several other recent publications.-"— 

The purpose of this work is to develop insight into how a detailed EOS model, and 
particularly ionization, affect spherically symmetric, targetless plasma liner implosions. We 
used the HELIOS-CR-i^ ID radiation-hydrodynamic code and the PROPACEOSi^ EOS and 
multi-group opacity database to repeat the simulation cases presented in Awe et ali^ In 
particular, we examined how the detailed EOS model affected the liner stagnation pressure 
Pstag (defined in Sec. Ill Cp . the stagnation time Tgtag (defined in Sec. Ill Cp . the implosion 
trajectory and minimum radius achieved, and scaling behavior of Pstag as a function of 
initial plasma liner parameters. In all cases, we used argon as did Awe et air- However, for 
a particular initial liner geometry, mass density, and velocity (corresponding to an initial 
energy of 375 kJ, relevant for PLX), we also ran simulations using H, D-T, ^He, ^Li, ^^B, Ne, 
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Kr, and Xe, in addition to Ar. In this work, as in Awe et al.— we set aside the important 
issues of 3D effects of the discrete merging plasma jets and 3D convergent instabihties 
such as Ray leigh- Taylor and others, which are addressed elsewhere,— and study spherically 
symmetric imploding plasma liners initiated at a jet merging radius Rm at which plasma 
jets are assumed to have merged. 

The remainder of the paper is organized as follows. Sec. [Ill describes our modeling ap- 
proach including the choice of liner initial conditions investigated, details of the codes used, 
and benchmarking against the results of Awe et ali^ Section IIIII presents the new simulation 
results including (1) comparison of results using local thermodynamic equilibrium (LTE) 
versus non-LTE EOS tables, (2) comparison of plasma liner implosions using a LTE tabular 
EOS versus a polytropic EOS, (3) Pstag scaling with initial liner density and velocity, and 
(4) different plasma species for one liner initial condition. Section |IV] provides a discussion 
and summary. 

II. MODELING APPROACH 

In this section we describe the simulation initial conditions, the codes used, and bench- 
marking of our own polytropic EOS results against the results of Awe et air- 

A. Simulation initial conditions 

The initial conditions used in this work are shown schematically in Fig. [T]and summarized 
in Table [H (columns 2-5). The rationale for these choices are described in detail in Awe et 
air- In short, they are intended to span a range of initial plasma liner velocities vq and kinetic 
energies KEq from PLX-relevant values (50 km/s, hundreds of kJ) at the lower end to what 
might be needed for an MIF standoff compression driver-"- (> 50 km/s, > tens of MJ) at 
the higher end. The initial liner ion density values tiq are what can be expected based on 
plasma jet densities that have already been achieved by the plasma guns designed and built 
by HyperV Technologies^ for use on PLX. In this work, we assume that discrete plasma 
jets have merged at Rm = 24.1 cm and formed a spherically symmetric imploding plasma 
liner with uniform mass density po? thickness 25.5 cm, temperature Tq = T^, = Ti = 1.0 eV, 
and uniform inward radial velocity Vq. The effects of non- uniform initial liner profiles are not 
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studied here but are a potential way of optimizing imploding plasma liner performance.— 
The void region (Fig. [1]) is modeled using very low density plasma with mass densities of 
10"^ kg/m^ (tabular EOS runs) and 10"^ kg/m^ (polytropic EOS runs), temperature of 
1.0 eV, and the same species as the liner. The different initial void mass densities led to 
slightly different computational zone setups but did not affect the overall liner evolution. 

Although this work does not focus on the physical evolution of a spherical plasma liner 
imploding on vacuum, which was described in detail in Awe et a/.,- we provide here a brief 
summary of its ID evolution: (1) the imploding liner converges toward the origin, and the 
liner density rises as R~^, (2) upon the liner reaching the origin, an outward going shock 
propagates radially outward, shock-heating the incoming liner material, and (3) the outgoing 
shock reaches the outer edge of the incoming liner, at which time the high pressure of the 
post-shocked liner material can no longer be sustained, and the entire system disassembles. 

B. Computational codes 

This work used the HELIOS-CR— ID radiation-hydrodynamic code with LTE and non- 
LTE EOS tables generated using the PROPACEOSi^ code. HELIOS-CR is a Lagrangian 
code that can use both a polytropic EOS (with adiabatic constant 7) or a tabular EOS. The 
HELIOS-CR code has been applied previously to the modeling of inertial confinement fusion 
(ICF) implosions.— PROPACEOS EOS and frequency- dependent opacities are computed 
using a combination of an isolated atom model and, at high densities at which cohesive 
effects are important, a quotidian-EOS-like model.— At the densities of interest in this 
paper, energies and pressures are computed for an isolated atom model in which atomic level 
populations are computed using Boltzmann statistics and the Saha equation (LTE case) or 
a collisional-radiative model (non-LTE case). The energies and pressures of the EOS tables 
are based on the free energy, thus providing for thermodynamically consistent energy and 
pressure values. For the non-LTE case, the effects of photoionization and photoexcitation are 
ignored when computing the temperature- and density- dependent tabular EOS data, as the 
radiation field is a non-local quantity, i.e., the radiation can originate in other parts of the 
plasma, and therefore is unknown when generating the tables. Detailed atomic models are 
utilized in which all ionization stages are included, e.g., for Ar, a total of > 10^ energy levels 
are included. More than 10^ transitions were used for the opacity calculations. Frequency- 
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dependent radiation transport is treated using a radiation diffusion model; fifty frequency 
groups were used in tliis work. Wlien modeling a plasma that is not optically thick (our case) 
with radiation diffusion, it is possible that radiation losses could be overestimated. Finally, 
we checked energy conservation at t = 5 i^is (near peak compression) for the polytropic, 
LTE tabular, and non-LTE tabular cases, using run 6 of Table |I] as an example. Energy 
is conserved to within 0.1%, 0.1%, and 1.1% for the polytropic, LTE, and non-LTE runs, 
respectively. These are very small effects compared to the differences we report in this paper, 
and thus we conclude that energy conservation does not impact our results. 

We previously verified^ a HELIOS-CR calculation against the convergent shock Noh 
probleroi^ and got better than 5% agreement for the highest resolution case we ran (50 ^m 
zone resolution over a sphere with 25 cm initial radius). We also did a grid resolution con- 
vergence study showing that an average grid resolution of 250 //m gave ~ 10% accuracy.- In 
this work, we used 400 computational zones over the initially 255 mm thick liner (average 
of 637.5 ^m/zone), which was chosen based on achieving a balance between getting reason- 
able agreement in the benchmarking studies described in Sec. Ill CI and the run time of each 
simulation (typically a few hours to a day per run). We used the automatic zoning feature 
of HELIOS-CR, and thus the zone size was much smaller than the average near boundaries 
and much larger elsewhere. 

When using a polytropic EOS, HELIOS-CR only allows use of a one-temperature (IT) 
model. The tabular EOS simulations allow for separate ion Tj and electron Te temperatures. 
All reported HELIOS-CR results with polytropic EOS are IT, and all reported results with 
tabular EOS are 2T with radiation diffusion. To verify that the use of IT versus 2T does 
not have a great impact on our results, we compared the time evolution of and Tj for the 
zone initially aX R = 25.1 cm and saw no difference. Furthermore, the Tj evolution from a 
IT run was virtually the same as that from the 2T run. 

C. Benchmarking against results of Awe et al. 

Because Awe et al.— used the ID Lagrangian radiation-hydrodynamics code RAVEN,— 
we first benchmarked HELIOS-CR against RAVEN by performing runs 1-8 of Table H] using 
a polytropic EOS and comparing with the RAVEN results. In this paper, all polytropic 
EOS runs used 7 = 1.6667. We also included thermal tranport and radiation diffusion to be 
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as consistent as possible with Awe et al.,— which found that in these spherically convergent 
simulations, non-physical temperature and pressure spikes occur near the origin due to 
compression of the leading edge material of the liner. Including thermal and radiation 
transport in the simulations helps prevent these non-physical spikes, which prevent the 
imploding plasma liner from reaching as small a radius as it otherwise would, leading to a 
much-reduced peak pressure.- 

We compared the quantities Pstag and Tstag, where Pstag is defined as the thermal pressure 
in the simulation zone that is initially 1.0 mm from the inner edge of the liner [R = 24.2 cm) 
averaged over a time duration T^tag, which is defined (and also shown graphically in Fig. |2J\) 
as the time duration between maximum pressure to the time when pressure sustainment 
ends due to the outward going shock reaching the outer edge of the incoming liner. We 
found that RAVEN P^tag results (Table [Tll) are 1.2-3.5 times higher than the HELIOS-CR 
results (Table HU]). The Tgtag results agree to within 20%. These differences were considered 
to be small enough to proceed with this study. 

We also compared HELIOS-CR and RAVEN simulation results for the time evolution 
(Fig. [2j\) of the pressure P, mass density p, and ion temperature Tj of a single computational 
zone (initially 1.0 cm from the inner edge of the liner, to be consistent with Fig. 4 of Awe et 
air-), and the radial profiles (Fig. [2B) of P, p, and electron temperature Te at a single time 
{t = 6.5 yus), for run 6 of Table HI Both comparisons indicate good qualitative agreement 
with some discrepancies. One difference is the increase in both P and Tj near the end of Tgtag 
for the HELIOS-CR but not the RAVEN results. The magnitude of the rise in P increases 
with increasing vq and is present throughout the entire stagnated liner, and is thus not a 
localized numerical effect in a single zone. The reason for this rise is not yet understood, but 
because this does not affect the conclusions of this paper, we have set the issue aside. There 
are also differences in the radial profiles at t = 6.5 /is, largely due to the slight difference in 
time evolution between the two simulations. At this time, there is an outward going shock 
for both the HELIOS-CR (at ~ 0.6 cm) and RAVEN (at ~ 0.5 cm) simulations, but the 
HELIOS-CR shock has propagated farther outward. The discrepancy in Tg inside the shock 
also appears to be due to the difference in time evolution, as slightly later radial profiles from 
the RAVEN results also fall toward the origin (see Fig. 4 of Awe et air-). These differences 
do not affect the conclusions of this paper. The rest of the paper focuses on comparing 
HELIOS-CR simulations using a polytropic EOS versus tabular EOS data. 
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III. RESULTS 



This section presents HELIOS-CR simulation results and analysis. Firstly, we discuss 
LTE versus non-LTE tabular EOS results, showing that there are only minor differences, 
thus justifying our choice to focus on the use of LTE for the subsequent results in the paper. 
Secondly, we compare argon plasma liner implosions between polytropic and a LTE tabular 
EOS to assess the effects of ionization and a more sophisticated treatment of EOS. Finally, 
we explore the effects of using different plasma liner species for run 6 of Table HI again using 
LTE tabular EOS data. 

A. LTE versus Non-LTE EOS tabular data 

The initial conditions we used (see Table H]) are in a parameter space in which LTE is 
not obviously correct. Non-LTE can have an effect on EOS {i.e., ionization/excitation) and 
radiation losses (excited states tend to be less populated). Below an argon ion density of 
~ 10^^ cm~^ at ~ 1 eV, the fractions of Ar I, Ar II, and Ar III as a function of temperature 
are noticeably different between LTE and non-LTE calculations. Above ~ 10^^ cm~^, cor- 
responding to later stages of liner implosion and stagnation, the EOS data are similar. We 
performed run 6 of Table [T] using both LTE and non-LTE tabular EOS data to explicitly 
assess the difference in results. Figure [3] shows a comparison of P, p, and Tj versus time 
for the computational zone that is initially 1.0 cm from the inside edge of the plasma liner. 
Note that for the tabular EOS runs, the pressure is calculated as njTj + UgTe, where and 
Tie are the ion and electron densities, respectively, Ue = Zrii, and Z is the mean charge state. 
The LTE and non-LTE results are virtually indistinguishable over the entire duration of the 
simulation, and there are also no appreciable differences in P, p, and T over the entire liner 
at several time steps (not shown). This indicates that either LTE or non-LTE EOS models 
could be used without altering the relevant simulation results for the purposes of this study, 
and thus we used LTE for all subsequent results presented in the paper. 

B. Tabular LTE EOS versus polytropic EOS 

To assess the effects of ionization and other effects included in the LTE tabular EOS 
data, we compared HELIOS-CR simulation results using LTE tabular EOS with those using 
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a polytropic (7 = 1.6667) EOS for runs 1-8 of Table HI For the computational zone that 
is initially 1.0 mm from the inside edge of the liner, the maximum pressure Pmax is 3.6-6.1 
times larger and Pgtag 3.9-8.6 times larger for the polytropic EOS results. Complete results 
for both the tabular LTE and polytropic EOS's are shown in Tables HI and IIIIl respectively. 
The Tj and p at peak compression were also higher in the polytropic EOS results as well 
(not shown). Figure H] shows P, p, Ti, and radial position versus time of the computational 
zone initially at 25.1 cm (1.0 cm from the inside edge of the liner) for the initial conditions 
of run 6 of Table [H The lower Pmax (~ zone thermal energy divided by zone volume) for 
the tabular EOS result is consistent with the combined effects of lower zone thermal energy 
(~ pT) and higher zone volume at peak compression, as compared to the polytropic result. 
Next, we investigate in more detail the reasons for the lower zone thermal energy and higher 
zone volume at peak compression for the tabular EOS run. 

Figure [5] shows various components of the total liner energy versus time for run 6 of Ta- 
ble [H The liner kinetic energy KE falls quickly around 5 ps when the leading edge of the liner 
reaches peak compression and stagnates. In the polytropic EOS case, the difference between 
the initial KE and the stagnated liner thermal energy TE is due mostly to radiative losses. 
In the tabular EOS case, around half of the initial liner KE goes to ionization/excitation 
energy, which is stored in free and excited electrons in the stagnated liner. The mean charge 
peaks around 7 (see Sec. IIII Dl) . The result of including the ionization/excitation physics is a 
drastic reduction in the liner compression at stagnation, leading to lower Pstag (see Tables [T] 
and [ml). 

C. Stagnation pressure scalings 

Awe et al— found the approximate scalings Pstag ~ ^0^''^ ^-iid Pstag ~ for the cases 
in Table [T] using a polytropic EOS model, and gave an heuristic argument for why Pstag ~ 
Vq ought to be an upper bound on the scaling dependence of Pstag on vq. The heuristic 
argument invoked the following assumptions: (1) perfect conversion of liner KE to liner TE 
at stagnation, (2) no radiation and thermal losses, and (3) no entropy production due to 
the outward going post-stagnation shock. Because the RAVEN runs of Awe et air- included 
radiation and thermal losses, the observed weaker Pstag ~ scaling is consistent with the 
heuristic upper bound of Pstag ~ v^. Our HELIOS-CR runs using a LTE tabular EOS violate 
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assumption (1) above due to the major ionization/excitation energy channel discussed in 
Sec. IIII B[ Thus, it is expected we should see an even weaker scaling than Pstag ~ ^'o^^"^- 

Figure [6] shows Pstag versus no and vq and fits to power law functions for HELIOS-CR 
simulations of the cases in Table I using a LTE tabular EOS model. While the fitted 
exponents are very similar across density and velocity groups in Awe et al.,— i.e., Pstag ~ 
^o.54±o.o2 Pg^ag ~ i)q'^^^°"®, our fitted exponents have a much wider spread, i.e., Pstag ~ 
^o.64±o.i4 Pgt^g ~ Vq^^^^'^'^. This larger spread in exponents is likely due to the fact 
that the LTE tabular EOS is sensitive to density and temperature in the imploding plasma 
liner, thus introducing a degree of variability that depends on the liner initial conditions 
not expected in the polytropic EOS simulations of Awe et al- Our Pstag ~ y^.giio.so gQ^jij^g 
(depending on no) is indeed weaker than the approximate Pstag ~ v^^'''^ scaling of Awe et 
al.,- as expected. Our Pstag ~ nQ^^^°'^^ scaling agrees with the Pstag ~ n^^"^ scaling of Awe 
et ali^ within our observed variation across velocity groups. 



D. Varying the liner species 

Due to the significant effect of ionization/excitation on achievable Pstag-, we explored the 
use of different plasma species for the same initial liner po and vq (and hence KE) of run 6 
of Table [B The only quantity that was varied was the initial number density no to keep po 
constant. We performed HELIOS-CR simulations with LTE tabular EOS data to investigate 
the effect of species on Pstag- The species chosen for study are those of interest as fusion 
fuel, blanket materials, or MIF implosion drivers for fusion energy applications: H, D-T 
(50%-50% by atom), ^He, ^Li, "B, Ne, Ar, Kr, and Xe. 

Figure [7JA. shows HELIOS-CR simulation results using a LTE tabular EOS model for 
P, volume, TE, and ionization/excitation energy at peak compression versus atomic mass 
number for the zone initially 1.0 cm from the inner edge of the liner. There is a clear 
trend of increasing peak P with atomic mass number, quantitatively consistent with the 
combined effects of decreasing volume and thermal pressure with atomic mass number at 
peak compression. The thermal and ionization/excitation energies are not monotonic with 
atomic mass number, and this is related to the different atomic and radiative properties of 
the different species. 

We also examined the mean charge and electron count versus time (Fig. [7)3) for the same 
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computational zone as above. Although the heavier species have more electrons to ionize, 
there are a larger number of the lighter species present by the ratio of their atomic mass 
numbers. It turns out that the peak mean charge for the heavier species reaches only about 
5-10 (Fig. [Tp), much less than the mass ratios between the heavier and lighter species, and 
thus the electron count is higher for the lighter species. Thus, the ionization/excitation 
energies are in general higher for the lighter species, in part explaining why they reach lower 
peak pressures for the initial conditions imposed. 

IV. DISCUSSION AND SUMMARY 

It appears that the key effect of using a LTE tabular EOS rather than a polytropic EOS 
is the significant diversion of liner KE into ionization/excitation, which is not included in 
the polytropic EOS model. The diversion of initial liner KE leads to less effective compres- 
sion, larger minimum volume at peak compression, and thus lower Pstag which is inversely 
proportional to volume. This is detrimental for reaching very high peak temperatures and 
pressures in targetless spherical plasma liner implosions, which was one of the design objec- 
tives of PLX. For a given initial liner KE, which is a natural experimental constraint, the 
results suggest the use of lower density and higher velocity to achieve higher peak pressure. 
The previous statement is confirmed by examining, e.g., the Pstag results from runs 3, 6, 
and 9 (375 kJ), or runs 4, 7, 10, and 13 (1.5 MJ) of Table [H Higher velocity should also 
be coupled with the use of lower atomic mass number species so as to fully strip the ions 
and eliminate the ionization/excitation and line radiation energy channels beyond a certain 
point in the implosion, after which the liner KE will predominantly go toward stagnation 
thermal pressure. This is in apparent contradiction with the results of Sec. IIIIDl which sug- 
gest the use of heavier species for higher Pstag but only because we constrained the species 
comparison study by enforcing the same Vq and po- In our species comparison study, the 
much higher uq of the lighter species leads to more of the initial liner KE going into ioniza- 
tion/excitation. Further study is needed for optimizing Pstag in targetless spherical plasma 
liner implosions for a given initial liner KE. 

The above discussion does not apply for the MIF application of a plasma liner compress- 
ing a magnetized D-T target plasma.- In that case, the imploding plasma liner is intended 
to be only a pusher that should complete its job of compressing the D-T target to fusion 
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conditions before the liner itself has a chance to get significantly compressed and heated. 
Thus, for initial liner no and vq that are likely to be limited by available plasma gun tech- 
nology, the MIF application likely will prefer a heavy liner species (to maximize pv"^) since 
the ionization/excitation should not matter until after the D-T target has reached peak 
compression. 

In summary, we have presented new ID HELIOS-CR radiation-hydrodynamic simulation 
results of spherical plasma liners imploding on vacuum using a LTE tabular EOS model. The 
results were compared in detail to HELIOS-CR simulations using a polytropic EOS model 
(7 = 1.6667), which in turn were benchmarked against the polytropic EOS simulation results 
of Awe et air- We found that the liner stagnation pressure Pstag is lower by a factor of 3.9-8.6 
using the LTE tabular EOS model. This is attributed to a significant amount of the initial 
liner kinetic energy going toward ionization/excitation of the plasma liner, resulting in less 
effective compression of the hner. We also studied the effect of different species on Pstag 
for the same initial liner geometry, mass density, and velocity. We found that with these 
constraints, the heaviest species reach the highest values of Pstag- This is attributed to less 
total energy going to ionization/excitation due to a smaller value at peak compression of 
the liner mean charge times the initial number density. 
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TABLE I. Initial conditions (same as Table U of Awe et ali^) and summary of HELIOS-CR argon 
simulation results for the indicated quantities, using LTE tabular EOS data. The results are for 
the computational zone that is initially 1.0 mm inside the inner edge of the liner. All cases have 
initial liner inner radius of 24.1 cm, liner thickness of 25.5 cm, and temperature of 1.0 eV. 
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Run Tstag (A's) Pstag (Pa) Pmax (Pa) PgiagTstag PstagTstag / K Eq 



1 


8.82 


2 37E+08 


3.17E+9 


2 OQE+S 


0.09 


2 


4.51 






1.40E+4 


0.15 


3 


2.29 


3 43E+lf) 


8 82E+n 


7.84E+4 


0.21 


4 


1.03 


6.29E+11 


4.26E+13 


6.48E+5 


0.43 


5 


8.84 


4.78E+08 


4.41E+09 


4.23E+3 


0.05 


6 


4.56 


6.79E+09 


1.32E+11 


3.09E+4 


0.08 


7 


2.33 


8.39E+10 


2.40E+12 


1.95E+5 


0.13 


8 


1.19 


8.32E+11 


2.35E+13 


9.86E+5 


0.16 



TABLE II. RAVEN argon simulation results using a polytropic EOS (7 = 1.6667) for runs 1-8 of 
Table U (reproduced from Table II of Awe et alX). The results are for the computational zone that 
is initially 1.0 mm inside the inner edge of the liner. 



Run Tstag (/"s) Pstag (Pa) 


Pmax (Pa) 


PstagTstag Pstag 


Tstag /KEq 


1 


10.30 


1.73E+08 


7.44E+08 


1.79E+3 


0.08 


2 


5.05 


1.87E+09 


1.09E+10 


9.44E+3 


0.10 


3 


2.55 


1.63E+10 


5.15E+10 


4.16E+4 


0.11 


4 


1.30 


1.78E+11 


3.19E+11 


2.32E+5 


0.15 


5 


11.15 


3.91E+08 


1.64E+09 


4.36E+3 


0.05 


6 


5.05 


4.11E+09 


2.71E+10 


2.08E+4 


0.06 


7 


2.65 


7.09E+10 


1.54E+11 


1.88E+5 


0.13 


8 


1.30 


3.86E+11 


7.70E+11 


5.02E+5 


0.08 



TABLE III. HELIOS-CR argon simulation results using a polytropic EOS (7 = 1.6667) for runs 1-8 
of Table m for benchmarking against the results of Awe et ali^ The results are for the computational 
zone that is initially 1.0 mm inside the inner edge of the liner. 
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FIG. 1. Schematic of the plasma liner initial conditions used in all simulations. The initial liner 
mass densities pQ and velocities vq are given in Table [H The initial liner temperature Tq = Tj = Tg 
is 1.0 eV. 
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FIG. 2. Benchmarking comparisons of (A) pressure P, mass density p, and ion temperature Ti (in 
the computational zone initially 1.0 cm from the inner edge of the liner) versus time; and (B) P, p, 
and electron temperature Tg versus radius at t = 6.5 ps between RAVEN and HELIOS-CR argon 
polytropic EOS simulation results of run 6 of Table [H 
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FIG. 3. Comparison of pressure P, mass density p, and ion temperature Tj (in the computational 
zone that is initially 1.0 cm inside the inner edge of the liner) versus time between LTE and 
non-LTE argon simulation results for run 6 of Table HI 
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FIG. 4. Comparison between argon LTE tabular EOS and polytropic EOS simulation results for 
run 6 of Table I for the time evolution of (A) pressure P and ion temperature Tj, and (B) radial 
position of the liner, for the computational zone that is 1.0 cm from the inside edge of the liner. 
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FIG. 5. Comparisons of the instantaneous total liner kinetic (KE), thermal (TE), ioniza- 
tion/excitation, and time-integrated radiated energies of run 6 of Table |T] using argon tabular 
LTE and polytropic EOS models. The latter does not include ionization/excitation. 
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FIG. 6. Pstag versus (A) tt-q and (B) Vq for the argon runs in Table U and corresponding power law 
fits. 
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FIG. 7. (A) Pressure, volume, thermal energy, and ionization/excitation energy, and (B) total 
electron count and mean charge versus liner atomic mass number at the time of peak compression, 
for the computational zone initially 1.0 cm from the inner edge of the liner, and for initial liner vq 
and po corresponding to run 6 of Table [H Liner species (from left to right on plot) include H, D-T, 
^He, ^Li, i^B, Ne, Ar, Kr, and Xe. 
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